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Abstract 

We calculate the nucleon parameters in nuclear matter using the QCD sum rules 
method. The radiative corrections to the leading operator product expansion terms are 
included, with the corrections of the order as beyond the logarithmic approximation taken 
into account. The density dependence of the influence of radiative corrections on the 
nucleon parameters is obtained. At saturation density the radiative corrections increase 
the values of vector and scalar self-energies by about 40 MeV, and 30 MeV correspondingly. 
The results appear to be stable with respect to possible variations of the value of Aqcd- 



1 Introduction 

In nuclear physics the nucleon parameters (self-energies) in nuclear medium are expressed in 
terms of meson exchange. In QCD sum rules (SR) approach to the problem the self-energies 
are expressed in terms of exchanges by systems of weakly correlated quarks. Until now these 
were just uncorrelated quarks. The correlations can be presented through the QCD radiative 
corrections. In the present paper we investigate the influence of these corrections on the values 
of nucleon self-energies in symmetric nuclear matter. 

The QCD sum rules (SR) method, suggested by Shifman et al. [Ij succeeded in expressing 
hadron characteristics in terms of expectation values of QCD operators. This approach was 
initially used for mesons. Later it was expanded by loffe et al. [2l E] to description of baryons 
(see also [Ij). In the SR method one considers the function n(g) which describes the propagation 
of the system with the quantum numbers of a hadron 



n(g) = I / ci4a;e^('?-)(0|Tj(a;)j(0)|0), 



with the local operator j{x) carrying the quantum numbers of the hadron. We shall consider 
the SR for proton. In this case 

n(g) = q^^^W{q') + m\q^) (2) 
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with 7^ and / standing for the Dirac and unit 4x4 matrices. 

The key point of the SR approach is the analysis of dispersion relations 



1 /• Im UUk^) , . 



(3) 



for the functions n*(g^). These equations are considered at large values of while < 0. 
The left-hand sides (LHS) of Eq. (3) are presented as power series of q~^, with the QCD 
condensates being the coefficients of the expansion. This presentation is know as the operator 
product expansion (OPE) [5J. The spectral densities Imn*(fc^) on the right-hand side (RHS) 
of Eq. (3) are usually approximated by the "pole-f-continuum" model 



Thus the position of the lowest pole m, its residue and the effective continuum threshold 
W'^ are the unknowns of the SR equations in vacuum. Usually the Borel transform of both 
sides of Eq. (3) is carried out. 

Several lowest order OPE terms have been found in [21 [3]. These were the contribution of free 
three-quark loop and the terms containing the vacuum condensates (0|gg|0), (Ol^G^^G^j.lO) 
and {0\qqqq\0), with q and G^^, the quark operators and the gluon field tensor. 

Later this approach was expanded for description of nucleons in nuclear matter O [71 [8] (see 
also [9] and references therein). In this case there are three structure of the polarization operator 



n„.(g,p) =2 / c?^xe^(^^)(M|j(x)j(0)|M) = g^7^n^(9,p)+PM7"n^(g,p) + /C(g,p), (5) 



with \M) and p - the vector of state of the matter and its four- momentum; p = (m, 0) in the 
rest frame of the matter. The dispersion relations for the functions ir{q,p) = Il^{q^, s) at fixed 
value of s = {p + qY 



are considered now instead of Eq. (3). The same approximations as in vacuum, i.e. the OPE 
and "pole+continuum" model were used for the LHS and RHS of the SR. The LHS contains 
now contributions of two types. There are expectation values of the same operators as in the 
case of vacuum, averaged over the ground state of nuclear matter {{M\qq\M) , etc.). There are 
also the QCD condensates which vanish in vacuum. In the leading OPE terms these are the 
vector condensate (M|g(0)7^g(0)|M) and the expectation value (M|g(0)7^D,yg(0)|M) caused 
by the nonlocal vector operator (M|g(0)7^g(x)|M). 

The in-medium characteristics of the nucleons, i.e. the vector self-energy J2v ^i-nd the Dirac 
effective mass m* are unknowns of the SR equations in nuclear matter. Two other unknowns 
are the residue at the nucleon pole and the continuum threshold, which obtain new values 
and W^. The SR approach [HI 13 HI [9] provided reasonable values of 200 MeV, 

m*—m ^ —300 MeV, which are consistent with the results obtained by nuclear physics methods. 
Similar results have been obtained in other modifications of the SR approach [TOl |TT] . 



Imm(g2) = \%5{k^ 



(4) 




(6) 
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Interactions between the quarks in polarization operators (1) and (5) manifest themselves 
in radiative corrections, which contain the powers of the QCD coupling constant as- In the 
asymptotics — oo the terms in which as is enhanced by the "large logarithm" In(g^) are 

the most important ones. The corrections (a^lng^)" have been obtained in all orders for the 
leading OPE terms [T^ and included into calculations carried out in [21 El El El E] • The SR are 
considered at finite values of the Borel mass M^, i.e. 

0.8 GeV^ < < 1.4 GeV^, (7) 

where is small enough for perturbative treatment of the radiative corrections. Perturba- 
tive calculation of such corrections in vacuum beyond the logarithmic approximation [T3] 
provided a numerically large coefficient of the lowest radiative correction to the leading OPE 
terms. Analysis of the role of radiative corrections in vacuum Borel transformed SR method 
which includes also corrections to the four-quark condensate [13] was carried out in p^. It was 
shown that at least for the standard current 

j{x) = eabc[ul{x)C'yf,Ub{x)]-f5-f^dc{x) 

these corrections modify mainly the value of the nucleon residue A^, while that of the nucleon 
mass suffers minor changes. In the present paper we are also using this current. 

Recent calculations [17] demonstrated that the lowest order radiative corrections to the 
leading OPE terms of the in-medium polarization operator (5) contain the coupling constant 
as multiplied by large coefficients, which are 7/2 in 11^ structure and 15/4 in 11^^ structure. 
Since ttg/Tr ?a 0.15 at the values of the Borel mass ^ 1 GeV^, new analysis of the sum rules 
with inclusion of radiative corrections is needed. 

We present the polarization operator 11^ determined by Eq. (5) as 

Um{q',s) = U{q^)+Up{q\s), (8) 

with the vacuum term 11 determined by Eq. (1) while Ilp{q^, s) describes interaction with the 
nuclear matter. We include radiative corrections to the vacuum term up to the contributions 
~ q~^ of OPE. We include also radiative corrections to the leading OPE terms provided by 
interaction with nuclear matter. We compare the results provided by inclusion of radiative 
corrections in several ways. 

It is well known that the sum of the terms (aglng^)" can be presented in terms of the 
function 

, Ino^/A^ , , 

) = to^w 

with A = Aqcd while fi = 500 MeV is the normalization point of the characteristic involved. The 
sum of corrections (a^lng^)" to a certain term multiplies it by a factor L~'^, with 7 reflecting 
the "anomalous dimension" of corresponding operator - see, i.e. \18\ . 

In the present paper we consider the Borel transformed sum rules for the nucleon in nuclear 
matter, considering the dispersion relations for the operator Up deflned by Eq. (8), employing 
the current j{x) suggested in |16j. We follow our paper |9j adding the radiative corrections to 
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the analysis carried out there. We include the logarithmic corrections and the correction of the 
order as beyond the leading logarithmic approximation (LLA) to the condensates of dimension 
d = 3, which are the vector and scalar condensates. We include also the LLA corrections to 
the main condensates of dimension d = 4, which are due to nonlocal structure of the vector 
condensate. 

We demonstrate that while the LLA corrections provide substantial changes of the values 
of vector and scalar nucleon self-energies and m* — m, these changes are to large extent 
compensated by the terms ~ a^, which are beyond the LLA. For example, at saturation value 
of nucleon density this composition of the radiative corrections adds 40 MeV and 30 MeV to S„ 
and m* —m correspondingly. These corrections effect mostly the values of the vector self-energy 
and of the nucleon residue. The values of m* — m and of the continuum threshold suffer minor 
changes. 

Our main numerical results are obtained for Aqcd = 150 MeV corresponding to ^^(l GeV^) = 
0.37 in one-loop approximation. We used these values in our previous papers on the subject. 
Somewhat larger values of Aqcd and GeV^) are often used nowadays - see, e.g. [1]. The 
possible modifications affect mostly the value of the nucleon residue, while the results for the 
self energies appear to be stable within several percent. We recall the results of inclusion of 
the radiative corrections to vacuum sum rules in Sect. II adding also some new data. The sum 
rules in nuclear matter with radiative corrections are considered in Sect. III. The results are 
discussed in Sect. IV. 



2 Radiative correction in vacuum 

Following [15] we present the lowest OPE terms of the operators 11'' (g^) and n^(g^) as 

^ = ^0 + ^4 + ^6 + ^; ^ = 53 + ^7 + ^9, (10) 

with lower indices showing the dimensions of the condensates, Aq is the contribution of the free 
3-quark loop. Perturbative inclusion of the lowest order corrections provides 



" 647r4^ fl^ \ 12 TT 2 TT fl^ 



2' 



^ 2 (0|gggg|0) / 5a, la, Q_ 
^ 3 g2 1^ 37r^/i2j' 

with = —q^ > 0. Radiative corrections to the terms and -By which are numerically small 
[21 [3], and to the terms A^ and Bg, which are known with poor accuracy are not included. 

The Borel transformed nucleon sum rules [21 [3] can be written as 

a{M\w^) =n''{M^), c\M\w^) = n\M^), (12) 
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with 

7^''(M2) = A^e-™'/*^', n\M^) = A^me-™'/*^' (13) 
describing the contributions of nucleon pole (A^ = 327r'^A^), while 

£^ = ^0 + ^4 + ^6 + ^, C = B; + Bj + B^. (14) 

Here Ai{Bi) denote the Borel transforms of the terms Ai{Bi) on the RHS of Eq. (10), multiplied 
also by 327r^. 

Explicit expression for Ai and Bi are given in |15]. Here we focus on the terms which 
contain the radiative corrections. Including the contributions (a^lng^)"' exactly and treating 
the corrections ~ beyond the logarithms perturbatively, we can write 

~ _ M^E2{W^/M'^)ro ~ _ 2a{M'^)M^Ei{W^/M^)rs 

ia - ^^re. (15) 

Here Ei{x) = 1 - (1 + x)e-^'; E2{x) = 1 - (1 + x + xV2)e~^. The function 

~, 2x InMVA^ . . 

comes from logarithmic corrections to the current j{x) and to condensates. Also in Eq. (15) 

a = -(27r)2(0|gg|0) ; = (27r)^(0|gggg|0) . (17) 

For i = 0,3 the factors 

ri = l^ Ci ; Co = — , C3 = - (18) 

TT iZ Z 

include the radiative corrections of the order beyond the LLA [H]. 

The situation with the contribution Aq is a bit more complicated. The four-quark condensate 
is presented usually under factorization assumption pQ 

{0\qqqq\0) = ((0|gg|0))^ a, = a'. (19) 

However, one should clarify the point, where Eq. (19) is postulated. Assuming, following [5] 
that the factorization a^lk^) = a?{k'^) to be true at certain k'^ = Mq, with Mq belonging to the 
interval determined by Eq.(7) and employing a(MQ) = a{iJ?')L'^/^ one finds in the LLA 

a,{Ml) = a'ifx')LiM',f'. (20) 

On the other hand dependence of the condensate 04 in framework of the factorization 
hypothesis was found in [H] 

a,{M') = a,{M',)(^l-^-^\n^Y (21) 



Of course, In ^ is not a "large logarithm" and hence is included perturbatively. Thus we can 
write 

Bs = 2a{^?)M'E,{W^/My,, A, = \a\^i^)L{Mlfl'r,{M\W^). (22) 



Here 



re(MM^^) = l-^(| + ln^ + ^(-^)), (23) 




with 

m = E 



includes the lowest correction beyond the LLA. 

In Eqs. (15) and (22) we assumed that the corrections of the order (a^lnM^)" and of the 
order compose independent factors. This can be justified by the fact that the former terms 
come from integration over momenta /i^ <^ /c^ <^ M^, while the latter come from fc^ ~ M^. 

In Table I we present the values of characteristics m, and W"^ obtained with the radiative 
corrections being included in various ways. We show the results with the radiative corrections 
totally neglected (i.e. all L = 1, = 1), the results in LLA based on Eq.(22) with L determined 
by Eq.(16), while = 1. We include the corrections beyond LLA, with determined by 
Eqs. (18) and (23). We show also the results for all perturbative inclusion of radiative corrections 
(PRIG)- see Eq.(ll). We assume the factorization point for the four-quark condensate (Eq.(19)) 
to be Mq = 1 GeV^. We present also the results corresponding to factorization point Mq = 
/i = 0.5 GeV obtained earlier in [15] . 



3 Radiative corrections in nuclear matter 

1. Sum rules without radiative corrections 

Following [HI El E] we consider the SR for 

n^ = A3p + A4p + A6p; nj = 53p + 56p; n^ = P3p + i^4p + P6p- (24) 

The Borel transformed SR are 

Cl{M\ Wl) = is, + A,, + is, = A„(M2) - A(M2); (25) 
£j(M2, Wl) = Bsp + Bep = m*A„(M2) - mA(M=^); (26) 
CliM^ Wl) = hp + P4p + Pep = -S,A„(M2), (27) 

with Am{M^) = A^e-™™/^^', A{M^) = A^e-^'/Af^ and are the effective value of the 
residue and of the continuum threshold in nuclear matter; can be presented in terms of 
and m* [9]. 
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The lowest dimension terms can be written as [9] 



A 



3p 



3p 



(pg)lnQ2; 

(M|Er7og*|M)g2ing2 

i 



(2^ 



{i = u, d) 
and 



■{N\Y.<tq'\N)p 



B 



3p 



(29) 



Since the SR are considered at s = + g)^ = const = 4m^ [SI [71 [H [9] , we must put {pq) 
(s - - g2)/2 in Eq. (28). From Eqs. (28), (29) we obtain [9] 

[s - m^)Eo{M^) - M^Ei{M^) _ 



As, = ~8n^M\J2q^^oq'\M)- 

i 

B^p = -A7C^N\J2q'q'\N)M'E,{M^)p; 

i 

Psp = ——{M\Y^q^^oq'\M)M'E,{M^ 



3m 



(30) 



The contributions of the fourth dimension A^p and come mainly from the nonlocality 
of the vector condensate {Q\Yl'f iS^)loD^q^{x = 0)|0) (there is also a small contribution Afp 

i 

caused by the in-medium gluon condensate - see [9j). Such contributions can be expressed in 
terms of the nucleon structure functions O El E] 

Aip = Alp + Alp ; 

1 

m Q 



^4p = ^ In — J da(^- aria{a) + 9rib{a)^p 



+ 



1 f Q''^ 



and 



In ^— I da 



{5a{pq') + 2a m )?7a(«) + 9m arib{a) p + 



(31) 



+ 



Gvr^m 



da In 



Q_ 



/2 



{-a{pq') - 2Q' )i]a{a) - 9m ar]b{a) 



(32) 



Here q' = q — pa, = — g^, Q'^ = —q'^, the functions 77a,fe(«) are defined following [H [9] 

{M\q{0)-fpq{x)\M) = 



^ J dae-'''^"'^\a{a)+ixpm J da e-'''^P''\b{a)\ P, 



with r]a{a) the standard deep inelastic nucleon structure function, the moments of the function 
(a) can be expressed in terms of those of the function r]a{(y) [H]. The cut off Ac will be 
eliminated by the Borel transform. 

For the Borel transforms of the left hand sides of Eqs. (31) and (32) (multiplied by 327r^) 
we find 

AlpiM"") = u'^iM^) + h^iM^y, P4p{M^) = uP{M^) + hP{M^), (33) 

with the contributions u"^ and coming from the first terms on the right hand sides of Eqs. 
(31) and (32), while h'^{M'^) and h^{M'^) come from the second terms. Expressions for u'^ and 
are presented in [9], while h'^ and require a special treatment. Due to the terms InQ'^/Q^ 
they have finite cuts, corresponding to the singularities in the u channel of the interaction 
between the baryon current and the quarks belonging to the nucleon of matter, i.e. by the 
exchange terms. In paper [H] these terms have been neglected. This was justified by claiming 
for the description of the nucleon in the Hartree approximation. 

Here we include the terms h'^ and h^. Performing the the Borel transform of 

, , (l+a)(g^ + A'^(a)) 

= 1? ■ 

with X^(a) = -^^{s—m^—rn^a) and using the numerical values of the moments of the structure 
functions rja [1.9J, we find the coefficients of expansion of the terms containing rja to be of 
the same order of magnitude. Thus the OPE in powers of series for these terms exhibits 
a poor convergence. Hence, following [8], we employ the explicit expressions 



hi(M^) = — f da ((s -m^ -2m'^a)Go(M\a) -GAM^a 
3m J ^ 


((-5(s - m^)a + 6mV)G'o(M2, a) + 
+ (4 + 5a)Gi(M^ a))r]a{a) - 18m^arib{a 



(34) 



k^iM^) ^ - Ida 



(35) 

Here Go{M^,a) = M2(1 -e-^'(")/^^'); Gi{M^,a) = (l - {1 + X\a)/M^)e-^'^"^/^'^') . We 
included only the lowest moments of the function r]h{a), which are related to the moments of 
the function ria{a) by the equations, presented in [8]. The values of these moments enable to 
expect the convergence of the latter expansion. 

The contributions A^p, B^p and Pep are determined by the 4-quark condensate. Here the 
calculations require certain model assumptions. We shall use the expressions presented in [9], 
based on the calculations of the four-quark condensates [20] in framework of the Perturbative 
Chiral Quark Model suggested in |2I]- Further development of this model is described in 
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2. Inclusion of radiative corrections 

Now we include the radiative corrections. If the LLA corrections (a^ In g^)*^ are included exactly, 
while the terms beyond the LLA are taken into account in the lowest order of perturbation 
theory, each term of Eq. (30), obtains a factor 

r = ^ (36) 

with X = q,I, X, and = 1 + ^c^. 

Following [ni [31 [7] we put ^ = Y = l'^ = 4/9; 7^ = 0. In the LLA 4^ = 0. Beyond the 
LLA we employ the recently obtained in- medium parameters c'' = 7/2 and = 15/4 [17] for 
corrections to the local vector condensate, while = 3/2 [Ti] . 

Considering the terms of the higher dimension, we include the factor L~^/^ corresponding to 
the anomalous dimension of the proton current. Note also that the nucleon structure functions 
employed in our paper reproduce their moments with the proper anomalous dimensions. 
Corrections to contributions of the nonlocal vector condensate beyond the logarithmic approx- 
imation require calculation of additional set of the Feynman diagrams. As it stands now such 
corrections are not known. 

We do not include the radiative corrections of the contributions A^p, B^p and P^p to the four- 
quark condensates, provided by the nucleons of the matter. This is because these contributions 
are obtained in framework of a certain model |2T], whose accuracy cannot be estimated. 

Thus we include the radiative corrections BLLA to the leading OPE contributions contain- 
ing condensates of dimension = 3. We include the LLA corrections to the terms containing 
condensates of dimension d = A. We do not include corrections to the terms with the con- 
densates of higher dimension. Such approach is possible since the leading OPE terms indeed 
provide numerically largest contributions. For example, considering Eq.(27) for and employ- 
ing the vacuum values of parameters and W^, presented in Table 1 we can see that at the 
saturation value of nucleon density po the local vector condensate contributes approximately 
330 MeV to S^, while nonlocal corrections and the four-quark condensate subtract about 50 
MeV and 70 MeV from this value. 

We compare several sets of results, corresponding to the cases, considered for vacuum in 
Sec.II. 

3. Numerical results 

Now we present the results of solutions of the SR equations. Following our previous works we 
put Q!s(l GeV^) = 0.37 which is consistent with the present data [23]. At the saturation density 
Po we find the vector self-energy = 230 MeV, while m* — m = —330 MeV, if the radiative 
corrections are totally neglected. 

Inclusion of the leading logarithmic corrections diminishes and m* — m by 70 MeV 
and 50 MeV correspondingly. Inclusion of radiative corrections beyond the logarithms makes 
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the values of and m* — m closer to those, obtained with total neglection of the radiative 
corrections. The data at saturation density are presented in Table II. Comparing the results 
of BLLA calculation and those obtained with perturbative inclusion of all ctg corrections, we 
see that they differ by about 20%. Thus it is important to include corrections (a^lng'^)" in all 
orders. Density dependence of nucleon parameters is shown in Figs. 1, 2. 

Somewhat larger values of ^^(l GeV^) are often used in nowadays calculations. The authors 
of employed as = 0.47 corresponding to Aqcd = 230 MeV in one-loop approximation, while 
the value as(lGeV^) = 0.55 corresponding to Aqcd = 280 MeV was used in [4J. We carried 
out calculations with these sets of parameters as well. Assuming ^^(l GeV^) = 0.55 we found 
very small changes of values of the vector self-energy S^, ^ 260 MeV and of the scalar one 
m* — m ^ —290 MeV. The residue suffered largest change, reaching the value 1.8 GeV^. 
Recall that in the case of vacuum sum rules variation of the value of also affected mostly 
the value of [15]. Dependence of the results on the value of Aqcd are shown in Table III. 

4 Summary 

We calculated nucleon self-energies and m* — m in nuclear matter in framework of the 
QCD sum rules approach, including the radiative corrections to the leading OPE terms. We 
demonstrated that there is a large compensation between the logarithmic corrections (a^ Ing^)" 
and the corrections of the order beyond the logarithmic approximation. We showed that at 
saturation value of nuclear density simultaneous inclusion of these corrections diminishes the 
values of the vector self-energy by about 40 MeV and of the scalar self- energy by about 
39 MeV. The threshold value suffers minor changes. The radiative corrections affect mostly the 
value of the nucleon residue A^. 

Result for the saturation density are presented in Table II. Density dependence of the 
parameters is shown in Figs. 1, 2. 

We analyzed also stability of the results with respect to variation of Aqcd and as{l GeV^) 
since other sets of values are often used in literature. We found the changes to affect mostly 
the value of the nucleon residue. The vector self energy S^, and the scalar self energy m* — m 
may vary only by percents at the saturation density. The results are given in Table III. 
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Table 1: Nucleon parameters in vacuum. Line 1-all radiative corrections are neglected. Line 
2 -radiative corrections are included in the leading logarithmic approximation (LLA). Line 3 
- corrections ~ ag are included beyond the LLA (BLLA). Lines 4 - results for perturbative 
inclusion of radiative corrections (PIRC 1) ag and ~ a^lng^ with factorization (19) assumed 
at = 1 GeV^. Line 5-(PIRC-2) with factorization assumed at Mq = 0.5 GeV 





m, GeV 


A2, GeV^ 


W^, GeV2 


Rad. corrections 








are neglected 


0.94 


1.84 


2.03 


LLA 


0.96 


2.13 


2.37 


BLLA 


0.91 


2.03 


1.96 


PIRC 1 


0.92 


1.94 


1.85 


PIRC 2 


0.95 


2.04 


1.91 



Table 2: Nucleon parameters at the phenomenological saturation value of nucleon density. 
Notations are the same as in Table. I 





S„, MeV 


m* — m, MeV 


Kn , GeV^ 




Rad. corrections 










are neglected 


229 


-329 


1.10 


1.72 


LLA 


160 


-380 


1.10 


1.77 


BLLA 


271 


-300 


1.41 


1.75 


PIRC 1 


334 


-251 


1.47 


1.80 


PIRC 2 


347 


-249 


1.48 


1.81 



Table 3: Dependence of the nucleon parameters on the values of Aqcd and ag at the phe- 
nomenological saturation value of nucleon density with radiative corrections included beyond 
the Leading logarithmic Approximation (BLLA). The results are presented for Aqcd = 0.23 
GeV, a,(lGeV2) = 0.47 and Aqcd = 0.23 GeV, ^.(IGeV^) = 0.47. The corresponding values 
of vacuum parameters are given in brackets. 



as{l GeV) 


E^, MeV 


m* — m, MeV 




w;i, Gcv" 


0.47 
0.55 


269 
264 


-291 (m = 0.93 GeV) 
-289 (m = 0.94 GeV) 


1.65 (2.35) 
1.81 (2.61) 


1.89(2.13) 
1.99 (2.26) 
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5 Figure captions 

Fig.l Density dependence of the vector self energy and of the scalar self energy m* —m. The 
nuclear matter density p is related to its saturation value po- Dotted lines: all radiation correc- 
tions are neglected. Dashed lines: radiative corrections are included in the Leading Logarithm 
Approximation (LLA). Solid lines: Corrections of the order as are included perturbatively 
beyond the LLA. 

Fig. 2 Density dependence of the nucleon residue and of the continuum threshold W^. 
Notations are the same as in Fig.2. 
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